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ABSTRACT 

We present the analysis of three Suzaku observations of a bright arc in the ROSAT All-Sky Survey 
1/4 keV maps at / w 247°, b « —64°. In particular, we have tested the hypothesis that the arc 
is the edge of a bubble blown by an extraplanar supernova. One pointing direction is near the 
brightest part of the arc, one is toward the interior of the hypothesized bubble, and one is toward the 
bubble exterior. We fit spectral models generated from 1-D hydrodynamical simulations of extraplanar 
supernova remnants (SNRs) to the spectra. The spectra and the size of the arc (radius « 5°) are 
reasonably well explained by a model in which the arc is the bright edge of a ~100,000-yr old SNR 
located ^1-2 kpc above the disk. The agreement between the model and the observations can be 
improved if the metallicity of the X-ray-emitting gas is ~l/3 solar, which is plausible, as the dust 
which sequesters some metals is unlikely to have been destroyed in the lifetime of the SNR. The width 
of the arc is larger than that predicted by our SNR model; this discrepancy is also seen with the 
Vela SNR, and may be due to the 1-D nature of our simulations. If the arc is indeed the edge of an 
extraplanar SNR, this work supports the idea that extraplanar supernovae contribute to the heating 
of the ^million-degree gas in the halo. 



Subject headings: Galaxy: halo 
background — 



— ISM: bubbles 
X-rays: ISM 



ISM: supernova remnants — X-rays: diffuse 



1. INTRODUCTION 



The discovery of shadows in the 1/4 keV soft X-ray 
background (SXRB) with ROSAT demonstrated that 
there is ~million-degree gas beyond the Galacti c disk , 
in the Galactic ha lo ([Burrows fc Mendenhalll 119911 : 
ISnowden et al.l Il991h . Subsequent a nalysis of data 
from the i?OS'^rAll-Sky Survey CRASS: IS^den et all 
Il997l) . and of spectra from pointed observations 
with XMM-Newton and Suzaku^ 
existence of this hot halo gas. 
gas is in coUisional ionization 
its t er nperature is ^1-3 x 10 ^ 
1998t iKuntz fc Snowde!^ [20001: 



has confirmed the 
Assuming that the 
equil ibrium (CIE) 
K (ISnowden et all 
lSmithetal.1 120071: 



Galeazzi et al.l 120071: iHenlev fc SheltonI 120081: iLei et all 
20091) . High-resolution X-ray absorption line spec- 
troscopy with the Cfeanrfra gratin gs also detected hot 
halo gas (e.g. lYao fc Wang! 120051 ). However, despite 
nearly 20 years of study, a fundamental question about 
the halo remains: how did the hot gas get there? 

Various mechanisms may be contributing to the hot 
halo gas. One possibility is that the hot gas originated in 
the disk, heated by stellar winds and supernovae (SNe), 
and w as transferred to t he ha l o via fountains or chimney s 
(e.g., iShapiro fc Fieldl [l976l : [N orman fc Ikeuchil [l989h . 
Another possible source of hot halo gas is gravitational 
heating of infalling intergalactic material, pr edicted by 
simulations of disk galaxy formation (e.g. I Toft et all 
I2002t iRasmussen et al.ll2009D . A third possibility, which 
is the subject of this study, is that the gas is heated 
in situ by SNe above the Galactic disk (Shelton 2006). 
X-ray spectroscopy is the key to distinguishing between 
these scenarios. For example, gas that has recently been 
heated by SNe will be underionized and gas that was 
heated by SNe in the distant past will be overionized 
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(e.g.. lSheltonl ll999^. while gas that has rapidly expanded 
out of the disk into the halo will be drastically ove- 
rioniz ed and recombining ([Breitschwerdt fc Schmutzled 
Il994[ ). Also, gas of extragalactic origin may have a dif- 
fe rent abundanc e pattern from gas of Galactic origin. 

ISheltonI ([2006D considered the soft X -ray emission from 
an ensemble of isolated supernova remnants (SNRs) of 
different ages and at different heights above the Galac- 
tic plane, taking into account the variation of the SN 
rate and ambient density as a function of height. She 
compared this emission to the 1 / 4 keV (R12) hal o emis- 
sion determined from the RASS. ISnowden et al.l ([1998D 
decomposed the observed 1/4 keV emission into a fore- 
ground component due to the Local Bubble (LB, a cavity 
in the interstellar medium of radius ^100 pc in which the 
Sun resides, thought to be filled with ~10® K gas), and a 
distant component, assumed to originate beyond the ma- 
jority of the Galaxy's H i. This distant component is a 
combination of halo emission and the extragalactic back- 
ground. Because the northern Galactic hemisphere con- 
tains such anomalous features as the North Polar Spur, 
and because at low latitudes one cannot clearly see the 
halo. lSheltonI ([2006D concentrated on high latitudes in the 
southern Galactic hemisphere. The average de-absorbed 
count-rate of the distant component for b < —65° is 
~800 X 10~^ R12 counts s~^ arcmin"^. Subtracting 
400 R12 co unts s~^ arcmin"^ for the extragalactic back- 
ground (e.g. lSnowden et al.lll998f ) l eaves ~4 0x 10 "^ R12 
counts s~^ arcmin"^ for the halo. ISheltonI ([2006D found 
that up to ~80% of this 1/4 keV flux could be explained 
by a population of extraplanar SNRs. In this scenario, 
the vast majority of the flux comes from old SNRs, cover- 
ing ~30-90% of the high-Galactic-latitude sky (the fact 
that they do not cover the whole sky explains the mottled 
appearance of the 1/ 4 keV halo emission in the maps of 
ISnowden et al.|[l998D . However, the individual old SNRs 
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a re rather dim, a nd so would be difficult to identify. 

ISheltonI ()2006D also calculated that ~1% of the high- 
latitude sky would be covered by young, bright remnants, 
which should be easier to identify. (Note that in this 
context, "young" means remnants that are still in the 
adiabatic phase - their ages could be up to ^10^ yr.) 
This fraction implies that ~1 young, bright extrapla- 
nar remnant is expected per Galactic hemisphere above 
|6| ^ 50°. SNRs have been found several hundred 
pc above the pla ne at low latitud es (e.g. , SN 1006 at 
z ~ 550 pc: .Winkler e t all l2003[), and IShelton et all 
()2007t l and iLei et all (|2009) have suggested that the X- 
ray and UV emission from an X-ray-bright region be- 
hind a nearby shadowing filament at I w 279°, b w —47° 
is consistent with the emission from a young remnant. 
However, to date, no isolated high-latitude extraplanar 
remnants have been confirmed. Here we examine an- 
other promising candidate for a wung extraplanar rem- 
nant, suggested by Shclton (20 Qa): a brig ht arc in the 
RASS 1/4 keV maps ( Snowden eiTall [19971 ) at I « 247°, 
b « —64°. This arc is the brightest arc- like feature in 
the southern Galactic hemisphere below ^—40°, and is 
shown in Figure [1] The shape of the arc is not due to 
absorption by intervening material. The right panel of 
Figure [T] shows the 1/4 keV data overlaid with contours 
showing the DIRBE -corrected IRAS 100-/im intensity 
(jSchlegel et al.l[l998h . which traces cool, absorbing mate- 
rial. One can see that there is increased 100- /im emission 
(and hence increased absorption of the background X- 
rays) to the lower left of the arc. However, the absorbing 
material does not completely follow the edge of the arc, 
implying that the arc's shape is not due to absorption. 

Although the RASS maps clearly show the presence of 
th e arc, ROSAT's low spectral resolution (E/AE ^ 1- 
3: ISnowden et al]|1997D makes detailed studies difficult. 
Therefore, we have used the X-ray Imaging Spectrometer 
(XIS; E/AE ~ 20 at £; ~ 1 keV: iKovama et al.ll2007i) 
onboard Suzaku (jMitsuda et al.ll2007D to obtain higher- 
resolution spect ra of the a rc and its surroundings, which 
we used to test ISheltonl 's ()2006[ ) suggestion that the arc 
is the edge of a young extraplanar SNR. Our three ob- 
serving directions are shown in the left panel of Figure [TJ 
One observation direction is near the brightest part of 
the arc ("On arc"), and one is toward the interior of the 
hypothesized SNR ("Arc interior"). The third observ- 
ing direction ("Arc background") is off the arc, outside 
the hypothesized SNR, and is intended to measure the 
ambient SXRB in the vicinity of the arc. Note that the 
XIS field of view (17.8' x 17.8') is much smaller than the 
radius of the arc (~5°), so we cannot obtain data from 
the entire arc region. 

The remainder of this paper is organized as follows. In 
^ we describe the Suzaku data reduction and spectral 
extraction. Our first goal is to establish if there are differ- 
ences in the halo component for our three observation di- 
rections. We investigate this question using CIE models, 
as described in fJS] We find that the halo is brighter and 
slightly cooler in the on-arc direction, compared with the 
other two directions. This CIE analysis provides a bench- 
mark against which to compare our subs eq uent a nalysis, 
using spectral models based on ISheltonl 's (|2006l ) hydro- 
dynamical simulations of extraplanar SNRs. This sub- 
sequent analysis, described in S]4l directly addresses the 
question of whether or not the arc is the edge of a extra- 



planar remnant. In addition, if the arc is the edge of a 
remnant, our analysis places constraints on the ambient 
density (corresponding to the height above the disk) and 
the age of the remnant. We discuss our results in [jS) In 
particular, we assess the SNR hypothesis on the basis of 
the spectrum, brightness, and gross morphology of the 
arc f ij5.2p . We finish with a summary in US) Throughout 
we quote Itr errors. 

2. SUZAKU DATA REDUCTION 
2.1. Initial Data Processing and Cleaning 

Table [1] shows the details of our Suzaku arc observa- 
tions. Our data were initially processed at NASA God- 
dard Space Flight Center (GSFC) using version 2 pro- 
cessing, specifically version 2.0.6.13 for the on-arc ob- 
servation, version 2.1.6.16 for the arc-interior observa- 
tion, and version 2.2.11.22 for the arc-background ob- 
servation. We carried out further data processing using 
HEAsoft^ version 6.6 and CIAO^ version 3.4, following 
guidelines available from the Suzaku Guest Observer Fa- 
cility at GSFC^. We used the set of calibration database 
(CALDB) files for the XIS released on 2009 February 
03, and the CALDB files for the X-ray Telescope (XRT) 
released on 2008 July 09. Throughout this paper we 
used only the data from the back-illuminated XISl chip, 
as it is more sensitive at lower energies than the front- 
illuminated chips. 

We first used the xispi tool to update the XIS gain 
calibration. We then cleaned and filtered the data. We 
selected events with grades 0, 2, 3, 4, and 6, and used 
cleansis to remove fiickering pixels. We excluded the 
times that Suzaku passed through the South Atlantic 
Anomaly (SAA), times up to 436 s after passage through 
the SAA, times when Suzaku^s line of sight was elevated 
less than 10° above the Earth's limb and/or was less 
than 20° from the bright-Earth terminator, and times 
when the cut-off rigidity (COR) was less than 8 GV. 
The thresholds used for the elevation of Suzaku^s line 
of sight above the Earth's limb and for the COR are 
higher than the defaults (which are 5° and 6 GV, respec- 
tively). The higher elevation threshold reduces contami- 
nation from the scattering of solar X-rays off the Earth's 
atmosphere, while the higher COR threshold reduces the 
particle background. We combined the data taken in the 
3x3 and 5x5 observation modes, and finally used the 
CIAO analyze_Ltcrv . si script to bin the 2.5-8.5 keV 
data into 256-s time bins and remove times whose count- 
rates differ from the mean by more than 3a. The result- 
ing cleaned XISl images are shown in Figure [2l 

2.2. Solar Wind Charge Exchange 

Observations of the diffuse soft X-ray emission from 
'^million-degree Galactic gas may be contaminated by 
geocoronal and heliospheric solar wind c harge exchange 
(SWCX) emission (Cravens 2000; Crave ns et aP [2001 
Robertson fc Cravens., 2003a .b: .Koutroumpa et al.ll2006l . 
20071) . Periods of enhanced SWCX emission have 



been associa ted with periods of increased solar wind 
proton flux (I Cravens et all 120011 : ISnowden et all 12004 
iFuiimoto et al.ll2007f) . 



^ http:/ /heasarc. gsfc.nasa.gov/lheasoft 
^ http://cxc.harvard.edu/ciao 

^ http: / /suzaku. gsfc.nasa.gov/docs/suzaku/analysis/abc/abc. html 
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Fig. 1. — ROSAT All-Sky Survey 1/4 keV maps of the southern Galactic hemisphere llSnowden et al.lll997l ). showing the arc that is the 
subject of this paper. Both panels show the same RASS data, which have been smoothed with a Gaussian whose standard deviation is 
2 times the pixel size. The units on the color bar are 10"'' counts s~^ arcmin"^. The coordinate grid shows Galactic coordinates. The 
small squares in the left panel show the Suzaku XIS field of view (17.8' X 17.8') for our three pointings (labeled "Arc interior", "On arc", 
and "Arc background"). The contours in the right panel indicate the DIRBE-corrected IRAS lOO-pm intensity llSchlegel et al.iri998l) . The 
contours are at 1, 2, 3, 4, and 5 MJy sr~^. 



TABLE 1 

Details of Our Suzaku Observations 



Observation 


Observation 
ID 


/ 

(deg) 


b 

(deg) 


Start time 
(UT) 


End time 
(UT) 


Usable exposure 
(ks) 


Arc interior 
On arc 

Arc background 


502070010 
502071010 
503104010 


253.29 
247.81 
240.49 


-62.74 
-64.51 
-66.01 


2008-01-15 19:09:14 

2007- 06-05 07:29:21 

2008- 12-30 06:07:54 


2008- 01-17 18:20:14 
2007-06-07 03:35:19 

2009- 01-04 08:19:23 


74.3 
72.6 
81.4 



Figure [3] compares the 0.3-2.0 keV XISl lightcurves for 
each of our three Suzaku observations with the contem- 
poraneous solar wind proton flux obtained from OMNI- 
Web^, which combines solar wind data from several dif- 
ferent satellites. In particular, the solar wind data for the 
arc-interior observation are from the Advanced Composi- 
tion Explorer {ACE), those for the on-arc observation are 
from ACE and Wind, and those for the arc-background 
observation are from Wind. 

In general the solar wind proton flux was fairly steady 
during our observations. The most notable exception 
to this statement is the first part of the arc-background 
observation, when the proton flux was greatly increased. 
However, there was no significant increase in the soft X- 
ray count-rate at this time. Despite the steadiness of the 
X-ray count-rate, we decided to err on the side of caution 
and removed times when the proton flux exceeded 2 x 
10^ cm~^ s~^, or when no solar wind data were available. 
These times are indicated by the gray datapoints in the 
lightcurves in Figure [3l The exposure times in Table [1] 
are those that remain after this additional filtering. 

http://omniwcb.gsfc.nasa.gov/ 



This procedure should help minimize contamination 
from bright, time-varying geocoronal SWCX. However, 
it should be noted that it is possible that heliospheric 
SWCX emission and some quiescent geocoronal SWCX 
emission still remain in our spectra. 

2.3. Point Source Removal 

We have previously found that automated source de- 
tection software d oes not work well on Suzaku images 
(Henley & Sheltori 2008), presumably because of the 
broad point spread function. We therefore identified in- 
dividual sources which may contaminate our diffuse spec- 
tra by eye from the 0.3-5.0 keV XISl images. Because 
the tool that we used to calculate non- X-ray background 
spectra (xisnxbgen version 2008-03-08) does not work 
correctly in sky coordinates, we worked in detector coor- 
dinates when identifying and removing individual sources 
from our XISl images. We used circles of radius 115 pix- 
els to exclude the events from around each identified 
source, except for one bright patch in the arc-background 
observation, for which we used a radius of 173 pixels. 
These radii correspond to w2' and «3' on the sky. The 
positions of the excluded sources are shown in the lower 
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Fig. 2. — Cleaned and smoothed 0.3-5.0 keV XISl images from our Suzaku arc observations. The upper row of images shows the data in 
sky coordinates. The lower row of images shows the data in detector coordinates. The black circles outline the regions that were excluded 
from the analysis - these circles were positioned by eye around sources that might contaminate the spectra of the diffuse emission. 

row of Figure O 

2.4. Spectral Extraction, Non-X-ray Background, and 
Response Files 

We extracted spectra from the full XISl field of view, 
excluding the aforementioned sources and times, and 
binned the spectra so there were at least 25 counts per 
bin. We used xisnxbgen to calculate non- X-ray back- 
ground spectra from the database of Suzaku night-Earth 
observations. For each observation we used the same 
parts of the XISl detector to extract the source spec- 
trum and to calculate the non- X-ray background spec- 
trum. We calculated redistribution matrix files (RMFs) 
using xisrmfgen, and ancillary response files (ARFs) 
using xissimarf gen, which takes into account the spa- 
tially varying contamination o n the optical blockin g fil- 
ters (OBFs) of the XIS sensor (l lshisaki et aIll2007D . For 
the ARF calculations we assumed a uniform source of 
radius 20'. 

3. SPECTRAL ANALYSIS; EQUILIBRIUM MODEL 

In this section we use a simple CIE model to look for 
variations in the halo emission in the vicinity of the arc. 
Our CIE model consists of unabsorbed emission from the 
LB and/or SWCX, absorbed thermal emission from the 
Galactic halo, and an absorbed power-law for the ex- 
tragalactic background due to unresolved active galactic 
nuclei (AGN). The model is described in more detail in 
H3.1\ and the results are presented in §3.21 Although this 
analysis does not directly address the question of whether 
or not the arc is the edge of an extraplanar SNR, it estab- 
lishes parts of the model that we use in our subsequent 
analysis using SNR models (iU, and it provides a bench- 
mark for comparison with SNR models. 



3.1. Model Description 
3.1.1. The Foreground Component 

When fitting a multicomponent model of the SXRB 
to a spectrum from a single direction, there can be a de- 
generacy between the foreground (LB/SWCX) and back- 
ground (halo 4 - extragalactic) e mission components. We 
therefore used iSnowden et aLl 's (|2000r ) catalog of shad- 
ows to fix the normalization of the foreground emis- 
sion in the vicinity of the arc. This catalog consists of 
367 shadows in the 1/4 keV RASS maps, five of which 
are in the vicinity of the arc: S2425M653, S2555M621, 
S2560M560, S2646M675, S2658M582. For each shadow, 
ISnowden et al.. (2000. ) decomposed the ROSAT Rl and 
R2 count-rates into foreground and halo count-rates 
(with associated errors). In the following test, we found 
that the foreground count-rates for three of these five 
shadows could usefully constrain the foreground emission 
model. 

We assumed that the foreground emission could 
be characterized with a IT CIE plasma model with 
T ~ 10^ K, as such a mode l provides a reason- 
able fit to the ROSAT dat a (jSnowden et all Il998l 
120001 : iKuntz fc SnowdenI 120001 ) . despite the fact that 
SWCX emission produces a different spectrum from 
thermal plasma emission . At first we fitted a IT 
Ravmond fc Smitlil (Il977l and updates) model with 
Anders fc Grevessg ()1989 ^ abundances to the foreground 



Rl a nd R2 count-rates for the five ,Snowdcn et al.l 
(|2000D shadows near the arc, we obtained a temper- 
ature of (S.S;?!) X 10^ K. This foreground tempera- 
ture is inconsistent with those found by previous stud- 
ies of the SXRB (T ^ 1.0-1.3 x 10^ K: ISnowden etall 



1998l[2000HKuntz fc SnoWdCTl[2000l : lGaleazzi et al.ll2007t 
Henlev fc SheltonI 120081) ." In addition, this foreground 
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5x10' 10= 1.5x10= 

Seconds since start of observation (2008 Jan 15 19:09:14 UT) 




5x10* 10^ 1.5x10^^ 

Seconds since start of observation (2007 Jun 05 07:29:21 UT) 




10= 2x10= 3x10= 4x10= 

Seconds since start of observation (2008 Dec 30 06:07:54 UT) 

Fig. 3. — XISl 0.3-2.0 keV lightcurvcs for each of our three 
observations, plotted alongside the contemporaneous solar wind 
proton flux. No solar wind flux data are available for the first 
~70 ks of the arc-background observation. The gray parts of the 
lightcurves correspond to times when solar wind data are missing, 



or when the proton flux exceeded 2 X 10 cm 
were removed from the data. 



These times 



model can be ruled out because it predicts an O viii Lya 



intensity of «9 photons cm ^ a 



(line units, L.U.), 



which far exceeds what is observed. A closer inspection 
of the shadows that we used showed that the foreground 
R2/R1 ratios of two (S2646M675 and S2658M582) are 
inconsistent with T ~ 1 x 10® K. When we eliminated 
these two shadows, we obtained a foreground emission 



model with T — 0.95 x 10® K and emission measure 
/ nidi — 0.0041 cm~® pc. We used this foreground 
model in our subsequent Suzaku analysis. 

3.1.2. The Halo and Extragalactic Components 

Although lYao fc Wan3 (|2007aD and iLei et all (|2009l ) 
have shown that isothermal and two-temperature halo 
models are inadequate for explaining all the available X- 
ray and UV emission and absorption data, these models 
still serve a useful purpose in characterizing the X-ray 
portion of the emission spectra. Our halo model con- 
sisted of a single CIE plasma component (IT model), 
which was adequate to characterize our data. The tem- 
perature and emission measure of the halo component 
were free parameters in the fitting. 

We modeled the extragalactic background as an ab- 
sorbed power-law . Its photon index was fixed at 1.46 
(|Chen et al.lll997[ ). but its normalization was a free pa- 
rameter. So that we could model the attenuation of the 
halo and extragalactic components, we obtained hydro- 
gen columns densities from th e Leiden-Argentine-B onn 
(LAB) Survey of Galactic H i (jKalberla et al.ll2005[ ) us- 
ing the HEAsoft nh tool. The hydrogen column densi- 
ties for our three observing directions are A'h — 1-9 x 
10^° cm~^ (arc interior), 1.6 x lO^*' cm^^ (on arc), and 
3.6 x 10^° cm~^ (arc background). 

3.1.3. Additional Details 

In order to better constrain the halo component at 
lower energies, we included Rl and R 2 (1/4 keV) data 
from the RASS (jSnowden et al.lll997t ). These data were 
extracted from circles of radius 1° centered on each of 
our Suzaku pointing directions using the sxrbg tool avail- 
able from HEASARC^. During the course of our analysis, 
we discovered a discrepancy between our Suzaku spectra 
and the ROSAT R45 (3/4 keV) count-rates - our models 
significantly underpredict the observed R45 count-rates. 
We will discuss this discrepancy in iJ5.1l We decided not 
to include the R45 data in our spectral analysis. We also 
did not include the R67 (1.5 keV) data, as this band is 
dominated by the extragalactic background, and it does 
not help constrain the Galactic therm al emission. 

We used XSPEC® version 11.3.2 (jArnaudl [1996^ to 
carry out the spectral analysis. For the thermal plasma 
components, we used the Astrop hysical Plasma E mis- 
sion Code (APEC) version 1.3.1 ("Smith et al ll2001i) for 
the Suzaku spectra and the [Raymond & S mithI ( 19771 
and upd ates) code for the ROSAT R12 data. We 
used the iRavmond fc SmithI code for the ROSAT R12 
data because APEC is inaccurate in that band, due 
to a lack of data on transitions from L-shell ions 
of Ne, Mg, A l, Si, S, Ar, and Ca.'^ The param- 
eters of each IRavmond fc SmithI component in the 
ROSAT model were tied to the parameters of the cor- 
responding APEC comp onent in the Suzaku model (see 
iHenlev fc SheltonI 120081 ). For the absorption, we used 
the XSPEC phabs model, which us e s cros s-sections from 
iBalucinska-Church fc McCam mon (1992"), with an up- 
dated He cr oss-section fromlY an et al. (1998). Through- 
out we used lAnders fc Grevessd (|l989D abundances. 

^ http:/ /heasarc.gsfc.nasa.gov/Tools/xraybg_help.html#command 
^ http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/xspecll 
^ http:/ /cxc. harvard.edu/atomdb/issues_caveats. html 
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We fitted the model to the 0.3-5.5 keV Suzaku + 
ROSAT R12 spectra, with the temperature and emis- 
sion measure of the halo components and the normaliza- 
tion of the extragalactic background as free parameters. 
The low-energy cut-off for the Suzaku spectra was chosen 
because the XISl calibration is uncertain below 0.3 keV. 
Although we did not expect much Galactic thermal emis- 
sion above ~1 keV, we included data up to 5.5 keV in 
order to constrain the extragalactic background. The 
high-energy cut-off was chosen to avoid the 5.9-keV Mn 
Ka line from the radioactive ^^Fe calibration source. 

3.2. Results for CIE Models 

We fitted our LB/SWCX + IT halo -t- extragalactic 
background model to each of our Suzaku + R12 spectra 
individually. The results are shown in Table [2] and in 
Figure HI Generally, the fits are reasonably good. The 
agreement between the model and the data is also good in 
the R12 band (not shown). The temperature of the halo 
is similar in all three directions ('^0.9-1.1 x 10^ K). How- 
ever, the halo emission measure is considerably larger in 
the on-arc direction - this is not surprising, given that 
the arc is brighter than its surroundings. 

The fact that the halo emission measure is larger in the 
on-arc direction than the arc-background direction sup- 
ports statement in the Introduction that the arc's shape 
is not due to absorption. If its shape were due to absorp- 
tion, we would expect the halo to have similar intrinsic 
intensities in these two directions - this is not what is 
observed. The halo is slightly cooler in the on-arc direc- 
tion compared with other two directions, although the 
difference in temperature is rather small. 

Although not relevant to the subsequent discussion, for 
completeness we note that the normalizations of the ex- 
tragalactic background in the three directions were 10.6 
(arc interior), 9.8 (on arc), and 8.1 (arc background) 
photons cm^^ s^^ sr^^ keV^ at 1 keV. 

4. SPECTRAL ANALYSIS: SUPERNOVA REMNANT 
MODELS 

The above-described analysis shows that there is in- 
trinsic variation in the halo's spectrum on and around 
the arc. However, by itself the above analysis has not 
helped us address the question of whether or not the arc 
is the edge of an extraplanar SNR. In this section, we 
test this hypothesis by fitting spectral models derived 
from h ydrodynamical simulations of extraplanar SNRs 
(jShelto n 2006) to our Suzaku -f R12 spectra. We look 
at SNR models with a range of ambient densities (cor- 
responding to different heights above the disk), ambient 
magnetic fields, and ages. We will assess these models 
on the basis of the spectrum, brightness and gross mor- 
phology of the arc. 

The hydrodynamical simulations from which our spec- 
tral models are derived are described in t j4.1l Our spec- 
tral model is described in i )4.2l In ii4.3l we present the 
results of the spectral fitting, and we identify the SNR 
models that are in best agreement with the observations. 
Finally, in i i4.4l we compare the radial surface brightness 
profile predicted by one of our best-fitting models with 
the observed ROSAT R12 and R45 profiles. 

4.1. Hydrodynamical Simulations of Extraplanar 
Supernova Remnants 
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Fig. 4. — Our observed arc-interior, on-arc, and arc-background 
Suzaku XISl spectra, plotted with the best-fitting CIE models ob- 
tained by fitting to the Suzaku + R12 data for each direction in- 
dividually. 



iSheltonI (|2006D carried out 1-D hydrodynamical simula- 
tions of extraplanar SNRs in 7 different ambient densities 
(no = 0.5, 0.2, 0.1, 0.05, 0.02, 0.01, 0.005 cm^^), corre- 
sponding to heights, z, above the midplane ranging from 
76 to 1800 PC (using the interstellar density model from 
iFerrierd I1998D . For each ambient density, she carried 
out simulations with (i?o/10^^ erg, Bos/mG) = (0.5,2.5) 
(which we call model type B), (0.5, 5.0) (model type C), 
and (1.0,5.0) (model type D), where Eq is the explosion 
energy and B^a is the effective magnetic field, which pro- 
duces a nonthermal pressure in addition to the ambient 
gas pressure. We have added a fourth type of model, with 
Eq = 0.5 X 10^1 erg and B^s = (model type A). The 
m odel paramete rs are summarized in Table[3](cf. Table 1 
in ISheltoiill2006( l. Table [3] gives the conversion between 
z and no, as well as model IDs that we will use below. 
The numerical part of the model ID indicates the height 
of the SNR in pc, while the letter (A-D) indicates Eq 
and i?cff. The models include thermal conduction, and 
the ionization evolution in the shocked gas is modeled 
self-consistently. 
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TABLE 2 
Spectral Fit Results 





Halo T 


Halo E.M.'' 






Observation 


(106 K) 


(10^3 cm-6 pc) 


xVdof 




Arc interior 


1.08 ±0.02 


17.8 ±0.6 


400.67/346 = 


1.16 


On arc 


0.95 ±0.02 


30.0t\i 


394.71/322 = 


1.23 


Arc background 


1.03 ±0.02 


23.3 ±i. 8 


371.31/339 = 


1.10 



^ Emission measure E.M. = J n^di 



TABLE 3 

SNR Model Parameters 



z 




Eq 


-Bcff 


Model 




(cm ^) 




(nCi) 


ID 


76 


0.5 


0.5 





76A 




0.5 


0.5 


2.5 


76B 




0.5 


0.5 


5.0 


76C 




0.5 


1.0 


5.0 


76D 


190 


0.2 


0.5 





190A 




0.2 


0.5 


2.5 


190B 




0.2 


0.5 


5.0 


190C 




0.2 


1.0 


5.0 


190D 


310 


0.1 


0.5 





310A 




0.1 


0.5 


2.5 


310B 




0.1 


0.5 


5.0 


310C 




0.1 


1.0 


5.0 


310D 


480 


0.05 


0.5 





480A 




0.05 


0.5 


2.5 


480B 




0.05 


0.5 


5.0 


480C 




0.05 


1.0 


5.0 


480D 


850 


0.02 


0.5 





850A 




0.02 


0.5 


2.5 


850B 




0.02 


0.5 


5.0 


850C 




0.02 


1.0 


5.0 


850D 


1300 


0.01 


0.5 





1300A 




0.01 


0.5 


2.5 


1300B 




0.01 


0.5 


5.0 


1300C 




0.01 


1.0 


5.0 


1300D 


1800 


0.005 


0.5 





1800A 




0.005 


0.5 


2.5 


1800B 




0.005 


0.5 


5.0 


1800C 




0.005 


1.0 


5.0 


1800D 



Note. — B^g = 2.5 and 5.0 /iG corre- 
spond to non-thermal pressures of 1800 and 
7200 cm~'^ K, respectively. 

Figure [S] shows how one such SNR evolves in time. 
T his mode l is mo del 1300B in Table [3l and remnant A 
m iSheltonI (fT999h : see that paper for more details. At 
early times (solid line, t — 25, 000 yr), the explosion cre- 
ates a hot cavity or bubble, bounded by the shock. The 
bubble produces copious soft X-rays - these X-rays come 
mainly from the hot dense region behind the shock, and 
the remnant is edge-brightened. As the bubble expands 
(dotted line, t = 100, 000 yr), it cools by adiabatic expan- 
sion. Nevertheless, the R12 emission remains bright and 
edge-brightened, although the R45 emission drops con- 
siderably. Between t = 100, 000 and 250,000 yr (dashed 
line), the dense rim of the remnant undergoes rapid ra- 
diative cooling, forming a cool, dense shell between the 
edge of the hot bubble at «100 pc and the shock at 
«105 pc. The X-ray count-rate drops considerably af- 
ter this occurs, although in this case the R12 emission 
remains edge-brightened. As time goes on, the hot bub- 
ble continues to cool (both adiabatically and radiatively) 
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Fig. 5. — The evolution of the density and temperature as a 
function of radius, and the predicted ROSAT R12 and R45 count- 
rates as a function of impact parameter, for SNR model 1300B (see 
text for details). 

and it gets progressively fainter in X-rays, while the cool 
shell widens. Eventually, the X-ray emission ceases to be 
edge-brightened (dot-dash line, t = 2 Myr), before the 
remnant eventually fades away altogether. 

Note that the times quoted above are specific to the 
chosen model - different model parameters (particularly 
the ambient density) will result in changes on different 
timescales. For example, a model with uq — 0.5 cm~'^ 
(corresponding to z = 76 pc) will form a cool shell be- 
tween t — 25,000 and 50,000 yr, and the R12 emission 
will cease to be edge-brightened hy t = 500, 000 yr. How- 
ever, the general pattern of change described above ap- 
plies to all the models. As described in the Introduction, 
Sheltonl (|2006f ) showed that the combined emission from 
an ensemble of such SNRs of different ages and at differ- 
ent heights could explain a large fraction of the observed 
1/4 keV halo emission at high Galactic latitudes. Here, 
we use the results of the SNR simulations summarized 
in Table [3] to calculate model spectra which we compare 
with our Suzaku + R12 spectra. 
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4.2. Spectral Model Description 

Each simulation dataset consists of hydrodynamical 
data and ion populations as a function of radius for up 
to 26 epochs, ranging from t = 2.5 x 10"^ to 2 x 10^ yr 
(the range of ages covered depends on the model). These 
data were used to calculate X-ray spectra. As the Suzaku 
field of view is much smaller than the size of the arc, 
we did not calculate spectra for the whole model rem- 
nant, but instead calculated projected spectra along var- 
ious sightlines through the rem nant. Fo r this p urpose, 
we used software developed bv ISheltoiJ (|1999( ). which 
uses the iRavmond fc SmithI ()1977( ) spectral code (up- 
dated by J. C. Raymond & B. W. Smith, 1993, pri- 
vate communication with R. J. Edgar) with abundances 
from lAnders fc Grevess^ (|1989| ) . The spectral calcula- 
tions take into account the (possibly non-equilibrium) 
ionization fractions output by the hydrodynamical sim- 
ulations. 

For each epoch of each SNR model, we calculated pro- 
jected spectra for sight lines through the SNR at various 
impact parameters (each was measured from the pro- 
jected center of the remnant). We normalized the model 
impact parameters such that the normalized value was 
for a sightline through the center and 1 for a sightline 
through the part of the rim at which the model R12 
count-rate was greatest. For example, for the model 
shown in Figure [S] the model impact parameters were 
normalized by dividing by 42.2 pc at t = 25, 000 yr, and 
by 72.7 pc at i = 100, 000 yr. We then used these sets of 
projected spectra as a function of this normalized impact 
parameter to create XSPEC table models®, one for each 
epoch of each SNR model. 

By eye, we estimated the center of the hypothesized 
bubble to be at (/,&) = (256.014°, -61.575°) (assuming 
that the bubble is a circle) , while the brightest part of the 
arc (in the R12 band) is at {l,b) = (247.90°, -65.09°). 
Therefore, the impact parameters of our on-arc and arc- 
interior sightlines, normalized to the impact parameter at 
which the R12 count-rate is greatest, are 0.936 and 0.341, 
respectively. We used projected SNR spectra calculated 
for these normalized impact parameters in our spectral 
analysis. 

As in ^J3l wc created a multicomponent emission model 
consisting of LB/SWCX, halo, and extragalactic compo- 
nents. However, instead of using a IT GIF model for the 
halo, we used a SNR -I- IT GIF model. The normaliza- 
tion of the brightness of the SNR component was a free 
parameter, i.e., the brightness of the model SNR emission 
could be scaled up or down to best match the observed 
spectrum. If a given SNR model accurately predicts the 
arc emission, its normalization should be 1. If the spec- 
trum predicted by the SNR component of our multicom- 
ponent model is too bright, the SNR model normaliza- 
tion will be less than 1. Similarly, if the model SNR 
spectrum is too faint, the SNR model normalization will 
be greater than 1. The normalization of the SNR compo- 
nent was constrained to be the same for the on-arc and 
arc-interior directions, and was set to zero for the arc- 
background direction (i.e., we assumed no SNR emission 
was contributing to the arc-background spectrum). The 
GIF portion of the halo component was constrained to 

* http: / / xspcc.gsfc.nasa.gov/docs / xanadu / xspcc / xspccl 1 / manual / : 



be the same in all three directions. As described below, 
the GIF component is needed to provide adequate fits to 
the data. This GIF component represents uniform hot 
halo gas that is not due a young SNR. This hot gas is of 
unknown origin, and could be due to old SNRs, outflows 
from the disk, or infalling extragalactic material. 

For the LB/SWGX emission, we used the model de- 
scribed in i)3.1.1l For the extragalactic background, we 
again used a power-law with a photon index of 1.46, the 
normalization of which was independent for each of the 
three observation directions. 

For each epoch of each SNR model, we fitted this 
LB/SWGX + GIF halo + SNR + extragalactic model to 
our three Suzaku + R12 spectra simultaneously. The free 
parameters during the fitting were the normalization of 
the SNR component, the temperature and emission mea- 
sure of the GIF halo component, and the normalizations 
of the extragalactic background. 

4.3. Results for SNR Models 

Here we evaluate our SNR models on the basis of 
their spectrum (which models give the lowest x^?), their 
brightness (which SNR components' normalizations are 
closest to 1?), and their gross morphology (which models 
predict an edge-brightened remnant whose radius is close 
to that of the arc?). 

Simultaneously fitting our LB/SWGX + GIF halo -|- 
SNR -I- extragalactic model to our three spectra yielded 
sets of best-fit parameters as a function of SNR age for 
each of the 28 SNR models. Figure [6] shows ^ 
function of SNR age for a subset of these SNR models. 
The results plotted are specifically for (a) model type A, 
with Eo = 0.5 X 10^^ erg and Beg = 0, and (b) model 
type B, with Eq = 0.5 x 10^^ erg and B^s = 2.5 ^G, 
but the other SNR models produce similar curves. The 
different curves correspond to SNe exploding at a series 
of increasing heights (that is, in a series of decreasing 
ambient densities). Note that not every model epoch 
yielded a valid fit - at later epochs the model SNR is 
brightest in the center, whereas the arc (if it is a SNR) is 
edge-brightened. As a result, not every epoch is shown 
in the curves in Figure [S] The GIE halo components 
from these fits have temperatures of '^0.8-1.1 x 10^ K 
and emission measures of ~0. 015-0. 03 cm~^ pc. The fits 
are considerably worse without this component. For ex- 
ample, without a halo component, the arc-background 
R12 count-rate is significantly underpredicted (602 x 
10~^ counts arcmin" , against an observed count- 
rate of (794 ± 13) X 10"*^ counts s^^ arcmin"^). With- 
out this component, the best-fitting model epochs give 
^ 1600 for 1012 degrees of freedom, against ^ 1200 
for 1010 degrees of freedom with the GIF component. 

For the lowest height {z = 76 pc) there are two minima 
in the curves, corresponding to SNR ages of 50,000 
and 150,000 yr. However, for z > 320 pc, there is a clear 
single minimum in the x^ curves, corresponding to SNR 
ages of 75,000-100,000 yr. (The exception is the model 
1800A, for which the best-fitting age is 250,000 yr.) 

The above-described curves tell us which SNR age 
gives the best fit to the spectra for a given SNR model. 
However, we would like to further discriminate among 
the models. To this end, we consider the normalizations 
idRjfetbftr^lNR components, the predicted radii at which the 
R12 emission is brightest, and x^- Figure [Tl^a) shows 
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Fig. 6. — Variation of with SNR age for SNR models at 
different heights, z. These results are from fitting the LB/ SWC X 
+ CIE halo + SNR + extragalactic model described in i|4.2l to 
our three Suzaku + R12 spectra, using (a) SNR model type A 
(Eo = 0.5 X 10^5^ erg, Beff = 0) and (b) SNR model type B {Eq = 
0.5 X 10^^ erg, B^g = 2.5 fiG). For clarity, the curves for z = 190, 
310, 1800 pc have been shifted upward by 50, 100,..., 300. 

the best-fitting normalization of the SNR component for 
each of the 28 SNR models. As was noted above, if a 
SNR spectral model accurately represents the arc emis- 
sion, then its normalization should be 1 - this is shown 
by a dashed line in the plot. Figure [Tl^b) shows the pre- 
dicted radii at which the R12 emission peaks. These radii 
were calculated for each SNR model using the best-fitting 
ages. The hydrodynamical simulations give these radii 
in parsecs. To convert these to angular radii, we calcu- 
lated the distances to the model SNRs using the nominal 
heights of the SNR models and assuming a Galactic lati- 
tude of —60°. If the center of the hypothesized bubble is 
at (/, b) = (256.014°, -61.575°) (see above), the observed 
radius at which the R12 emission peaks is w5° - this is 
shown by a dashed line in the plot. Figure [ZIc) shows 
the best-fit for each model. 

For z — 76 pc, the best-fitting SNR models are ~6-12 
times too bright, and a factor of ~4 too large. For z = 
190 and 310 pc, the predicted radii are within a factor 
of 2 of the observed value, but the best-fitting models 
are ^20-120 times too bright. The exception is model 
190D, which is less than 3 times too bright. However, its 
predicted radius is 2.5 times too large. For z — 480 pc 
the predicted radii are in very good agreement with the 
observed value, but the best-fitting models are still an 
order of magnitude too bright. 

The best agreement between the models and the obser- 
vations is for z > 850 pc. In terms of the normalization of 
the SNR component, the models with Eq — 0.5 x 10^^ erg 
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Fig. 7. — The best-fitting (a) normalization of the SNR compo- 
nent, (b) radius of brightest R12 emission, and (c) for each of 
the 28 SNR models that we investigated. The results are plotted as 
a function of the height of the original SN explosion, corresponding 
to different ambient densities (see Table 1 in Shclton 2006). The 
letters A-D denote models with different explosion energies, Eq, 
and effective ambient magnetic fields, B^ti (see Table [Sjl. The hor- 
izontal dashed lines denote a SNR normalization of 1, and a radius 
of brightest R12 emission of 5° (see text for details). Note that the 
values of x^ are for one fewer degree of freedom than in Figure |6] 
because for each model we have selected the epoch with the lowest 



(model types A, B, and C) are better than those with 
Eo = 1.0 X 10^1 erg (model type D). The best-fitting 
Eq — 0.5 X 10^^ erg models are less than a factor of 
4.5 too bright, and the radii are within a factor of 2.2 
of the observed value (the exception is model 1800A, 
which is a factor of ~10 too faint). Apart from this 
one model, all models for z > 850 pc give a best- 
fitting SNR age of 100,000 yr. However, we cannot eas- 
ily discriminate among the models for z > 850 pc with 
Eq = 0.5 X 10^^ erg, apart from ruling out model 1800A. 
For example, models 1800B and 1800C have the best- 
fit SNR normalizations closest to 1, but these models 
also have larger values of than the models for 850 or 
1300 pc. Also, at a given height, all the models have 
similar values of x^- 

Figure [8] shows our three Suzaku spectra along with 
the best-fitting spectral model obtained from the SNR 
model 1300B {Eq = 0.5 x lO^^ erg, B^s = 2.5 /iG). This 
model is in the middle of the range of acceptable heights, 
and in the middle of the range of magnetic fields used. 
The model generally does a reasonable job of fitting the 
spectra (including the R12 data, which are not shown), 
although the ~0.6-0.9 keV flux is underpredictcd in all 
three spectra (in the arc-background spectrum, the fiux 
is also underpredictcd at lower energies) . Over the whole 
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Fig. 8. — The three Suzaku spectra with the best-fitting 
LB/SWCX -I- CIE halo + SNR + extragalactic model obtained 
using SNR model 1300B (Eq = 0.5 X 10^1 erg, B^g = 2.5 ^lG). 
The individual model components are also plotted. 

Suzaku band, the SNR component is brighter in the on 
arc direction. However, above 0.6 keV the SNR com- 
ponent is brighter in the arc-interior direction, implying 
that the X-ray emission is harder toward the center of 
the SNR than toward the edge. 

4.4. Comparison of the SNR Model with R12 and R45 
Profiles across the Arc 

An additional comparison we can make between our 
model and the observations is to look at profiles of the 
ROSAT R12 and R45 count-rate across the arc. Such 
a comparison is made in Figure [9l which shows (a) R12 
and (b) R45 profiles along a great circle on the celestial 
sphere starting at the estimated center of the hypothe- 
sized SNR at I = 256.014°, b = -61.575°, and passing 
over the arc, as near as possible to the three Suzaku ob- 
servations. The model profiles were calculated using the 
best-fit parameters obtained in the previous section using 
SNR model 1300B {Eq ^ 0.5 x lO^^ erg, B^s = 2.5 ^iG). 
For each direction a long the model profil e, we used the 
column density from lKalberla et "all (|2005D calculated us- 



Distance from estimated center of SNR (deg) 

Fig. 9. — Profiles of the observed (a) R12 and (b) R45 count-rates 
across the arc (crosses) compared with the profile predicted by our 
best-fitting LB/SWCX + CIE halo -|- SNR -I- extragalactic model 
obtained using SNR model 1300B {Eq = 0.5 X 10^1 erg, B^g = 
2.5 fiG) (solid line). The individual model components are also 
plotted. The horizontal axis shows the distance from the estimated 
center of the hypothesized SNR at I = 256.014°, b = -61.575°. 
The arrows show the positions of the three Suzaku observations 
(left to right: arc interior, on arc, arc background). 

ing the HEAsoft nh tool. Although we found some varia- 
tion in the normalization of the extragalactic background 
among our three Suzaku spectra (see §3.2p . we do not 
know how the extragalactic background varies in the di- 
rections between our Suzaku observations. Therefore, for 
the purposes of Figure [9] we assumed a constant value of 
9 photons cm^^ s^^ sr^^ keV~^ at 1 keV. The conclu- 
sions drawn from Figure [9] are unlikely to be affected by 
relaxing that assumption. 

Figure (H^a) shows that the variation in the model 
R12 intensity across the arc is due both to the edge- 
brightened SNR, and to variations in the CIE halo in- 
tensity (which are due to variations in iVu, as this model 
component does not vary intrinsically) . The model over- 
predicts the R12 intensity near the center of the hypoth- 
esized bubble, and underpredicts the width of the arc. 

In the R45 band (Figure ^h] ) , the model underpre- 
dicts the observed intensity by about a third (~30- 
40 X 10~^ counts arcmin"^). For the middle data 
point, the discrepancy is even larger, but this data point 
may be contaminated by two inaccurately removed point 
sources (see i )5.1[ below). The disagreement between the 
model and the observations shown in Figure [^b) is not 
due to a problem with the SNR model itself, but instead 
is indicative of the general discrepancy between Suzaku 
and the ROSAT R45 band aheady noted in gSTSl 

5. DISCUSSION 

The main goal of this project is to test the hypothesis 
that the bright arc in the 1/4 keV SXRB at I w 247°, 
b « —64° is the edge of a bubble blown by an extraplanar 
SN. In H5.2I below, we will discuss this SNR scenario. 
However, first we discuss the discrepancy between the 
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TABLE 4 

Comparing Model and Observed R45 
Count-rates 







Observed 


Direction 


Model 


0.5" 


1" 


Arc interior 


75 


109 ±9 


119 ± 5 


On arc 


69 


168 ± 12 


122 ± 5 


Arc background 


55 


100 ±9 


95 ± 5 



Note. — All values are in 10~^ R45 counts 
arcmin"^. The observed count-rates were 
averaged over circles of radius 0.5° and 1°. 



Suzaku and ROSAT R45 intensities, noted in §fJ3TTT2]and 

HH 

5.1. The Discrepancy between Suzaku and the ROSAT 
R4-5 Band 

As was noted in H3.1.31 during the course of our analy- 
sis we discovered a discrepancy between our Suzaku spec- 
tra and the ROSAT R45 (3/4 keV) count-rates. Ta- 
ble m compares the R45 count-rates predicted by our 
Suzaku fit results with the observed count-rates, aver- 
aged over circles of radius 0.5° and 1°. For the R45 
count-rates extracted from the 1° circles, the models un- 
derpredict the rates by ^40 x 10~^ counts s~^ arcmin~^ 
for the arc-interior and arc-background directions, and 
by 53 X 10~^ counts s~^ arcmin^^ for the on-arc di- 
rection. The on-arc discrepancy increases to 99 x 
10"^ counts s~^ arcmin"^ when we use a 0.5° circle. The 
SNR model also underpredicts the R45 count-rate ( §4.4p . 
The fact that the CIE and SNR models both exhibit this 
discrepancy implies that the discrepancy is not due to a 
p roblem with the SNR mo del. 

iHenlev fc SheltoiJ (|2008l ) also noticed a discrepancy 
between their Suzaku spectra and the corresponding 
ROSAT count-rates. They partially overcame this dis- 
crepancy by adding a vphabs absorption component to 
their model. This extra component modeled contamina- 
tion on the XISl optical blocking filter over and above 
the contaminatio n alre ady included in the CALDB. 
IHenlev fc SheltonI (|2008f l found that they needed an ex- 
tra carbon column density of 0.28 x 10^* cm~^, in ad- 
dition to the CALDB value of 3.1 x 10^^ cm'^ at the 
center of the XISl chip. As the systematic uncertainty 
on the contamination thickness is ^^0.5 x 10^^ cm~^,^ this 
correction is not unreasonable. However, we were unable 
to obtain a good fit simultaneously to the Suzaku, R12, 
and R45 data for the on-arc direction, even with an extra 
vphabs absorption component in our model. 

We think that the discrepancy may be partly due to 
the fact that the R45 emission is mottled on scales of 
^^1° (see Figure [T0[b]). The observed R45 count-rates 
are averaged over bright and faint mottled regions. For 
the on-arc direction, this includes a particularly bright 
region, which may be partly due to inaccurate removal 
of the point sources IRXS J023800. 5-390505 and IRXS 
J023734. 5— 391925 (whose positions are shown in Fig- 
ure [T0|) . Because of the smoothing in Figure [TOl and 
because the size of the XIS field of view is similar to 
the RASS pixel size (17.8' x 17.8' versus 12' x 12'), one 

^ http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/xisO. html 



cannot accurately determine the R45 count-rate in the 
area exactly corresponding to the XIS field of view. If 
our Suzaku pointings happen to be toward fainter parts 
of the mottling, while the ROSAT count-rates include 
both bright and faint parts of the mottling, then the 
ROSAT fluxes will be systematically brighter than the 
Suzaku fluxes. In contrast, although the R12 data is 
much more variable over the whole sky than the R45 
data. Figure [Toy a) shows that the R12 emission does not 
seem to be as mottled on small angular scales. 

The discrepancy may also be partially due to quies- 
cent SWCX emission that is at a higher level in the 
ROSAT data than in the Suzaku data. Variations in 
the SXRB count-rate on a timescale of a few days, re- 
ferred to as "long-term enhancem ents" (LTEs), wer e 
removed from the RASS data (Snowden et al.l [1995?). 
These LTEs are now thought to be due to variations 
in the heliospheric a nd/or geocoronal SWCX emission 
([Cravens et al.l [2001[ ). However, even after the LTEs 
have been removed, a quiescent level of SWCX emis- 
sion may remain in the data. In the R45 band, this 
SWCX emission would be dominated by O vii and O viii 
emission. At high ecliptic latitudes, such as that of 
the arc (/? ^ —60°), the heliospheric SWCX emission 
from these lines is expected to be brighter at solar max- 
imum (such as 1990/1991, when the RASS was car- 
ried out) than at solar minim um (when our Suzaku ob - 
servations were carried g ut) (|Koutroumpa et all 120061 ). 
iKoutroumpa et all ()2007[ ) give O Vli and O VIII inten- 
sities for various Chandra, XMM-Newton, and Suzaku 
observations of the SXRB carried out at solar maxi- 
mum and solar minimum. Using the intensities predicted 
by their "ground level" model (that is, excluding short- 
term solar wind enhancements) for the so-called "south- 
ern Galactic fllament" (SGF), we estimate that the R45 
count-rate at high ecliptic latitudes due to quiescent he- 
liospheric SWCX is ~10 X 10~^ counts s~^ arcmin"^ 
higher at solar maximum than at solar minimum. This 
could explain ~l/3 of the discrepancy between Suzaku 
a nd R45. 

lYoshino et all ()2009f ) compared observed ROSAT R45 
count-rates with those predicted by Suzaku spectra from 
14 different directions (1 direction, the North Ecliptic 
Pole, was observed twice). They found that the ROSAT 
intensities were systematically brighter than the Suzaku 
intensities. For 5 of these directions, the ROSAT count- 
rates may be signiflcantly conta minated by LTEs. Af - 
ter removing these 5 directions, lYoshino et al.l ()2009l ) 
found that the ROSAT rates are an average of 17 x 
10~^ counts s~^ arcmin"^ higher than the correspond- 
ing Suzaku rates. They found that much of this offset 
could be due to a difference in the point source sensi- 
tivity between the two datasets. They also considered 
variations in the heliospheric SWCX emission between 
the two datasets. 

Figure [TT] shows a comparison of the observed ROSAT 
R45 count-rates and the r ates predicted by Suzaku; the 
plot shows the results from lYoshino et all ((2009, see their 
Figure 6) and our arc results. For the purposes of this 
plot, we re-extracted the observed R45 count-rat es using 
a circle of radius 0.3°, to match lYoshino et al.l ([2009). 
As can be seen, the discrepancy that we see for our 
three spectra i s not unusually large when compared with 
lYoshino et al.l 's results. 
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Fig. 10. — Maps of the RASS (a) R12 and (b) R45 count-rate in the vicinity of the arc. The data have been smoothed with a Gaussian 
whose standard deviation is 2 times the pixel size. The units on the grayscale bar are 10~® counts arcmin"^. The squares shows 
the approximate XIS fields of view for our Suzaku pointings (top to bottom: arc interior, on arc, arc background). The surrounding 
circles are of radius 0.5° and 1°. The two smallest circles near the on-arc pointing show the sources IRXS J023800. 5— 390505 and IRXS 
J023734. 5— 391925, mentioned in the text. The white line outlines the outer edge of the arc. 
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Fig. 11. — Comparison of the R45 count-rates predicted by var- 
ious Suzaku observations with the observed ROSAT R45 count- 
rates (after Yoshino ct al. 2009). The crosses show data from 
lYos hino ct al.|. The d atapoints surrounded by squares are identi- 
fied by Yos hino et al.l as being from possibly contaminated regions 
of the RASS. Our arc data are shown by open circles. The diagonal 
solid line shows equality. 

Our foregroun d LB/SWCX model is derived from 
RASS R12 data (jSnowden et al.ll2000t ). and so our fore- 
ground model should already include a reasonable esti- 
mate of the contamination of the R12 emission by solar 
maximum heliospheric SWCX emission and unresolved 
point sources. This model does not make a major con- 
tribution in the Suzaku band (see Figs. |4] and [8]). In 
addition, as noted above, the R12 data do not seem to 
be mottled on small angular scales as the R45 data (see 
Figure fTO]) . We therefore think that our combining the 
Suzaku and R12 data should yield reasonably accurate 
results. 

5.2. Is the Arc the Edge of a Supernova Remnant? 

In we used a multicomponent model of the SXRB, 
which include emission from a SNR, to analyze our three 
Suzaku -I- R12 spectra. In i i4.3l we showed that our spec- 
tra and the observed size of the arc are reasonably well 
explained by a model in which the arc is the bright edge 



of a -100,000-yr old SNR, blown by a SN explosion with 
Eq = 0.5 X 10^^ erg at a height of ~l-2 kpc. (Despite its 
age, we consider such a SNR to be young, as it is still in 
the adiabatic stage of its evolution, before the formation 
of a cool shell.) 

It is important to note that the SNR model described 
in 12] does not fit the spectra better than the GIF models 
described in [JH As the three fits in Table [2] are com- 
pletely independent, we can add the values of and the 
degrees of freedom to give = 1167 for 1007 degrees of 
freedom for the GIF model. In constrast, our preferred 
SNR models give ~ 1220-1250 for 1009 degrees of 
freedom. (Note that our two models are not "nested", 
so we cannot use the i^-test to measure the significance 
of this difference in x^-) The SNR models appear to do 
slightly worse in terms of x^, but we point out that the 
SNR scenario provides a good explanation for the arc's 
general morphology. 

It should also be noted that there are some discrepan- 
cies between the SNR models and the observations: the 
SNR models are generally either a factor of '-^3-4 too 
bright or a factor of ^2 too small. In addition, the mod- 
els do not match the observed width of the arc (i i4.4p . 
We will first discuss some possible explanations for these 
discrepancies fi^ ij5.2.1l and l5.2.2)) . before discussing other 
possible observations of the arc ( §5.2. 3|) . 

5.2.1. The Brightness and Size of the Arc 

At early epochs, when the model SNRs are still edge- 
brightened, their X-ray surface brightnesses and physical 
sizes are most strongly affected by the ambient density, 
riQ. As the X-ray emission comes from shock- heated am- 
bient medium, a denser ambient medium gives brighter 
remnants, as well as hastening the formation of the dense 
shell. A denser ambient medium also leads to smaller 
remnants, as in the Sedov phase the radius at a given 
age is proportional to n^^^^ (e.g. ISpitzerlll978l §12. 2b). 
The apparent size of an object also depends on its dis- 
tance. In the case of our SNR models, the distance is 
derived fron i the height, z , corresponding to the ambi- 
ent density (jSheltonl 120061 ) . using a Galactic latitude of 
-60°. 

Model 1800G (Eq = 0.5 x lO^i erg, B^g = 5.0 /xG, 
no = 0.005 cm~^) gives the best result in terms of the 
SNR surface brightness - of all our SNR models, the 
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best-fitting normalization for 1800C is the closest to 1 
(0.78t[!;S]^). However, this model gives a SNR radius that 
is a factor of ~2 too small. A simple way of resolving 
this discrepancy is to assume that the SNR is a factor of 
~2 closer than expected, without changing the ambient 
density. The surface brightness and physical size of the 
SNR would be unaffected, while the apparent size of the 
SNR would be ^2 times greater, and in better agreement 
with the observed size of the arc. This suggested resolu- 
tion requires the height at which ng = 0.005 cm"^ to be 
«1 kpc, as opposed to 1.8 kpc. Is a density of 0.005 cm"^ 
at z — 1 kpc plausible? 

Th e density model no{z) used bv lSheltonI (|2006l ) comes 
from iFerrierd (|1998f ). At z = 1 kpc, this model gives 
no(H -I- He) = 0.016 cm"^^ assuming 10% He by num- 
ber. This density is dominated by the neutral medium 
(37%) and the warm ionized medium (61%), t he remain - 
ing 2% being due to the hot ionized medium. iFerrieref s 
model for the warm ionized mediu m comes from the dis - 
tribution of free electrons derived bv' Cordes et al.l (|1991h , 
while her model for the neutral medium comes from 
iDickev fc LockmanI ()1990D . 

We used a Monte Carlo method to investigate whether 
or not a density of 0.005 cm"^ at z = 1 kpc was plau- 
sible. For each trial, we randomly varied each param- 
eter of the warm-ionized and neutral models according 
to its error bar, and calc ulated no(z = 1 kpc ) from the 
resulting density model. iDickev fc LockmanI ()1990f ) do 
not quote errors for the neutral model's parameters, so 
we assumed 50% errors. We included the hot ionized 
medium in no(z), but did not vary its model parameters. 
We carried out 100,000 trials. Approximately 24% of the 
trials were discarded because the random number gener- 
ator gave one or more negative (and hence unphysical) 
model parameters for those trials. Of the remaining tri- 
als, 7% gave a density no{z = 1 kpc) < 0.005 cm"^. We 
therefore cannot rule out the possibility that the density 
at z = 1 kpc is as low as 0.005 cm"^^ corresponding to 
model 1800C. As a result, we cannot rule out the possi- 
bility that the arc is well-described by model 1800C at 
z ~ 1 kpc instead of z = 1.8 kpc, although this model is 
a less favored option. 

The models with Eq = 0.5 x 10^^ erg for z = 850 pc 
(850A, 850B, 850C) are in better agreement with the 
observed radius of the arc, but they are a factor of 
~3 too bright. As stated above, the surface bright- 
ness of a model SNR depends on the density of the 
ambient medium. However, as the emission is domi- 
nated by line emission from ionized metals, the SNR sur- 
face brightness also depends upon the metallicity of that 
medium. If the metallicity of the interstellar medium 
above the di sk were a factor of ^3 lo wer than our as- 
sumed value (|Anders fc Grevessil 19891) , the model SNR 
surface brightness would be in better agreement with 
that of the arc. We also note that uncertainties in the 
model line emissivities may be important, although this 
is more difficult to quantify. 

The halo gas-phase abundances of Si, Mg, and Fe, 
which are all important line emi tters in the 1/4 keV band , 
are ^ 0.6, ~0.3, and ^0.2 solar (jSavage fc Sembachlll996l . 
using lAnders fc Grevessel[l989l as a reference). However, 
these values do not imply a subsolar metallicity for the 
halo; instead, these elements are assumed to be depleted 
on to dust. In contrast, S has a solar gas-phase abun- 



dance in the halo (jSavage fc Sembachlll996f) . implying it 
is not depleted on to dust. Also, the abundance of Ne 
(which, being a noble gas, is not expected to be depleted) 
toward an X- ray binary (4U 182 0—303) at z « 1 kpc is 
1.2 ±0.2 solar (|Yao fc Wang||2006f ) . (However, this source 
is at low Galactic latitude [b = —7.91°], so the sightline 
samples the disk as well as the halo.) 

These measurements suggest that the halo has a solar 
metallicity, at least for z < 2 kpc. However, the X-ray 
emission from a halo SNR depends upon the gas-phase 
metallicity in the shock-heated gas. Whether or not this 
metallicity is subsolar depends on how quickly the dust 
is destroyed behind the shock. In gas of temperature 
T = 1 X 10® K and density n, the rate of chang e of a 
dust grain's radius a due to thermal sputtering is (jSeabl 
fT987i) 

and so the lifetime r of a dust grain in the gas is 
da/dt 

The density n = 0.1 cm ^ used in the above expres- 
sion is the approximate density in the immediate post- 
shock region of the SNR models at z = 850 pc; this 
region is where most of the X-ray emission originates. 
The lifetime of a 100-A dust grain in such a SNR is an 
order of magnitude greater than the SNR age given by 
our spectral analysis (^10^ yr), and dust grains larger 
than 100 A would survive even longer. We would there- 
fore expect that the X-ray emission from a young halo 
S NR would reflect th e depl eted gas-phase abundances 
in iSavage fc SembachI (|1996[ ) , not the total halo abun- 
dances. 

The effect that these depleted abundances would have 
on our model SNR spectra depends on the details of 
which elements are depleted and by how much. However, 
as the gas-p hase halo abundances of s everal elements are 
~l/3 solar (jSavage fc SembachI [l996l) . it seems reason- 
able t o suggest that our model SNR spectra, calculated 
using lAnders fc Grevessi (|1989f ) abundances, may be a 
factor of ~3 too bright. Using a metallicity of ~ 1/3 so- 
lar for a SNR with Eq = 0.5 x lO^^ erg at z = 850 pc 
(riQ = 0.02 cm~3) would bring the model SNR surface 
brightness into better agreement with that of the arc. 
If the arc is a young extraplanar SNR whose X-ray- 
emitting gas has a metallicity of ~l/3 solar, our analysis 
would then favor a height z ~ 1 kpc for the arc, as a SNR 
at greater height (lower ambient density) would have too 
low a surface brightness, as well as being too small. 

5.2.2. The Morphology of the Arc 

As noted in i i4.41 the observed radial width of the arc 
is larger than that predicted by our SNR model. We in- 
vestigated this discrepancy by looking for known Galac- 
tic SNRs at an analogous stage in their evolution (i.e. 
roughly half way to the formation of the cool shell), in 
order to see how their observed X-ray morphologies com- 
pared with our model. 

We carried out a SNR simulation with Eq = 0.5 x 
gj.g^ j^^^ _ Q (]-nodel type A), with an ambient den- 
sity (no — 2 cmT^) representative of the disk, rather 
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than the halo. From this simulation, we estimated that 
a ^15,000-yr old SNR in the disk would be at a similar 
stage in its evolution as our best-fitting halo SNR model 
for the arc. We found that the Vela SNR (G263.9-3.3) 
provides a good analog to the arc: based on the ap- 
parent origin of X-ray-emitting explosion fragments be- 
yond the blastwave and the proper motion of the Vela 
pulsar (PSR B0833-45). lAschenbach et al.l (|1995f ) esti- 
mated the age of the remnant to be 18, 000 ± 9, 000 yr, 
consistent wit h the pulsar's spin-down age of 11,000 yr 
(|Tavlor et a l.| [r993l . 

As with the analogous models of halo SNRs, the new 
simulation with uq = 2 cm~^ predicts that a SNR of the 
age of Vela would have an X-ray-bright rim whose width 
is -1/10 of the radius of the SNR (cf. Figs. [5] and [H). In 
contrast, Vela exhibits an asymmetrical bright rim to the 
north and east whose wi dth is —half the remnant radius 
(|Aschenbach et al.|[l995l ). 

Our SNR models underpredict the widths of the arc 
and of the Vela SNR's X-ray-bright rim. These discrep- 
ancies may be due to the limitations of our SNR sim- 
ulations, which are one-dimensional, and assume a uni- 
form ambient density and a uniform non-thermal pres- 
sure in the radial direction. Three-dimensional simula- 
tions could include non-uniform ambient densities and 
different magnetic field geometries, which could affect 
the predicted X-ray morphology. Such simulations could 
also simulate hydrodynamical instabilities, which would 
broaden the apparent width of the arc. E. Raley (2009, 
private communication) has provided us with the results 
of a 3-D simulation of a SNR with Eq = 0.5 x 10^^ erg 
evolving in zero m agnetic field at z = 400 pc (see also 
iRaley et al.ll2007l ). Using these results, K. Kwak (2009, 
private communication) has provided us with spectra cal- 
culated for various sightlines across the remnant, assum- 
ing CIE. For a young remnant, the 3-D model does indeed 
predict a broader X-ray-bright rim than the 1-D model. 
However, the 1-D SNR models that we used include self- 
consistent modeling of the ionization evolution, which we 
used to calculate non-equilibrium X-ray spectra. Devel- 
oping a 3-D SNR model that includes ionization evolu- 
tion is beyond the scope of this paper. 

We conclude this discussion of the arc's morphology 
by noting that the arc does not trace a full circle. Such 
asymmetries are not uncommon in SNRs (see, e.g., the 
collection of SNR images at the ROSAT Guest Observer 
Facility^*' and the Chandra Supernova Remnant Cata- 
log^^). The arc's asymmetry may be due to a non- 
uniform ambient medium and/or a complicated magnetic 
field geometry, both of which are beyond the reach of our 
1-D simulations. 

5.2.3. Non-X-ray Observations, and Future X-ray 
Observations 

We have considered further ways in which we could test 
the hypothesis that the arc is the edge of an extraplanar 
SNR. Our X-ray spectral analysis implies that, if the arc 
is a SNR, then it is young, in the sense that it has not 
yet formed a dense cool shell. This is unfortunate, as it 
means we do not expect Ha emission due to recombina- 
tion in the cooling shell, nor H i emission from the gas 

http:/ /heasarc. nasa.gov/docs/rosat/gallery/snrs. html 
http:/ /hea- www. cfa.harvard.edu/ChandraSNR/ 



that has already cooled. The Souther n Ha Sky Survey 
Atlas (SHASSA; iGaustad et an[20?)l iFinkbeineii [20031 ^ 
shows no sign of an Ha-emitting shell around the arc. 

The hydrodynamical SNR models that we have used 
in our spectral analysis also make predictions for the 
column densities of various highly ionized metals (e.g., 
N V, O VI, O VII, O viii). In principle, therefore, far- 
ultraviolet and X-ray absorption line spectroscopy could 
be used to test the results of our spectral analysis. Un- 
fortunately, in practice it would be extremely difficult to 
detect unambiguously an enhanced ion column density 
due to an extraplanar SNR, given the sparsity of stars in 
this region. 

Table shows various ion column densities predicted 
for a 100,000-yr old SNR at z > 850 pc with Eq = 0.5 x 
10^^ erg. The predicted column density is dependent 
upon z and BcS- models at lower z (and hence higher 
ambient density) or with larger B(.g give larger column 
densities. Table [5] shows the range of values predicted by 
the various models. The N v and O vi columns (and, to 
a lesser extent, the O vii column) are peaked at the edge 
of the remnant (see Figs. 8 and 9 in lSheltoiilll998t that 
model corresponds to our model 1300B). For these ions. 
Table O shows both the peak column, and the column in 
the SNR interior. The O viii column does not peak at 
the edge of the SNR. 

Table [5] also shows the range of column densities ob- 
served in multiple directions across the sky. In general, 
the column densities predicted by the SNR are smaller 
than the spread in the observed column densities, and 
so an increased column density in the vicinity of the arc 
could not be unambiguously attributed to an extraplanar 
SNR. Furthermore, the predicted O vii columns are of 
the same order of magnitude as the errors on the O vii 
measurements. Some of the models (with z = 850 pc) 
predict peak N v and O vi column densities near the rem- 
nant edge that are significantly larger than the variation 
across the sky. However, observing these large column 
densities would require a fortuitously positioned back- 
ground source. In general, using absorption line spec- 
troscopy to test whether or not the arc is the edge of an 
extraplanar SNR would require column density measure- 
ments from several different sightlines on and around the 
arc. 

Perhaps the best way of further testing the hypothesis 
that the arc is the edge of an extraplanar SNR will be 
with future, higher-resolution X-ray spectrometers. Al- 
though our Suzaku spectra are consistent with CIE mod- 
els, it is possible that higher resolution spectra could 
reveal unambiguous signs of non-equilibrium ionization 
that we would expect behind the blastwave of a young 
extraplanar SNR. Sensitivity in the 1/4 keV band, where 
the arc is bright, but where we currently only have the 
low-resolution ROSAT R12 data, would be particularly 
useful. 

5.3. Other extraplanar SNRs? 

IShelton et all (|2007f ) and iLei et all (|2009( l have ana- 
lyzed the halo O VI and X-ray emission in the direction 
I « 279°, b « —47°, which samples an X-ray-bright re- 
gion visible in the 1/4-keV RASS maps. They used a 
nearby {d = 230 pc) shadowing filament to separate the 
foreground (LB/SWCX) and background (halo) emis- 
sion. This filament is visible toward the upper-right cor- 
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TABLE 5 

Predicted SNR Ion Column Densities and Observed Ion Column Densities 



Predicted SNR Column Density Observations 

Peak Interior Column Density IN umber ot 

Ion (lO^-* cm-2) (1014 cm-2) (ioi4 cm-^) Sightlines Ref. 



N V 0.2-1.6 0.03-0.19 

O VI 3.9-14.2 0.6-2.1 

O VII 13-43 10-23 

O VIII ■ ■ ■ 1.3-8.5 



0.57 ±0.36=* 9 1 

2.5±1.0'' 91 2 

~100-200 3 3,4,5 

~10-100 2 3,4 



References. — ( 1) ISembach fc Savagj 11992: ('21 ISavage et al.l 12003: (3) 
lYao fc Wanil2007al : f4) IYao fc Wanall2007bl : (5) lYao et aL l [200l 

Note. — The predicted column densities are from 100,000-yr old SNRs at 2 > 
850 pc with Eq = 0.5 x 10^^ erg. The observed values indicate the variation in the 
column densities detected from multiple sightlines over the sky. 
^ Mean it standard deviation 



n er of Figure [J 

IShelton et all (|2007f) showed that the O vi-to-1/4 keV 
emission ratio was a good age diagnostic for an extra- 
planar SNR. For the halo beyond the shadowing fila- 
ment, they found that this ratio was consistent with that 
from a ~40,000-70,000-yr old SNR with no = 0.01 cm'^, 
Eo = 0.5 X 10^1 erg, B^ff = 2.5 (model 1300B in Ta- 
ble[3]). More generally, they found that if the hot gas that 
they observed can be compared to that in an undisturbed 
extraplanar SNR, the time since heating is ~10*-10^ yr. 

Using a combi nation of Suzaku , ROSAT R12, and 
FUSE O VI data, iLei et all (|2009D constructed a differ- 
ential emission measure (DEM) model for the halo be- 
yond the shadowing filament. Their DEM is a broken 
power-law between T ~ 10^ K and T ~ lO'^ K, with a 
break at Tbrcak 10^ K. They suggest that the lower- 
temperature part of the DEM (T < Tbreak) is due to an 
extraplanar SNR. The O vi and soft X-ray intensities of 
this lower-temperature component are best matched by a 
-180,000-yr old SNR with uq ^ 0.02 cig-^ (z = 850 pc). 
Given the uncertaint i es in the modeling, IShelton et al.l 
(|2007f ) and iLei et al.l (|2009t ) obtained consistent results: 
both studies found that the O vi and 1 /4 keV X-ray emis- 
sion from the halo beyond the shadowing filament are 
consistent with the emission from a pre-shell-formation 
S NR at z ~ 1 kp c. 

ISheltonI (j2006t ) predicted that there would be approx- 
imately one bright, pre-shell-formation SNR per hemi- 
sphere at high-latitudes. Our arc analysis and the anal- 
ysis o f the filament region (|Shelton et al.|[2007t iLei et ah! 
l2009l ) suggest the presence of two halo SNRs at similar 
heights and of similar ages in the southern Galactic hemi- 
sphere. Given the uncertainties in the computer simula- 
tions and SN r ate, this res ult sh ould not be considered 
a refutation of ISheltonl 's (|2006( l prediction. However, 
given this result, we expect that most other directions 
in the southern Galactic hemisphere would not exhibit 
the properties of young, extraplanar SNRs. This could 
be tested by measuring the O vi-to-1/4 keV ratio for a 
large number of high-latitude directions. 

6. SUMMARY 

We have analyzed a set of three Suzaku spectra ob- 
tained from pointings on and around a bright arc that ap- 
pears in the 1/4 keV soft X-ray background at I « 247° , 
b « —64°, and have tested the hypothesis that the arc 
is the edge of an extraplanar SNR. For this purpose we 



have used spectral models generated from 1-D hydrody- 
namical simulations of SNRs at a variety of heights above 
the disk. We also supplemented our Suzaku spectra at 
lower energies with ROSAT R12 (1/4 keV) data. 

The three Suzaku + R12 spectra and the observed size 
of the arc are reasonably well explained by a model in 
which the arc is the bright edge of a '--^100,000-yr old SNR 
at a height z ^ 1-2 kpc (ambient density no = 0.02- 
0.005 cm^'^), blown by a supernova with Eq = 0.5 x 
10^^ erg. The remnant is still in the adiabatic stage of 
its evolution, having not yet formed a cool shell. 

The agreement between the model and the observa- 
tions is not perfect: the models are generally either a 
few times too bright, or a factor of ~2 too small. The 
agreement between the model and the observations can 
be improved if the metallicity of the X-ray-emitting gas 
is ~l/3 solar. Such a subsolar metallicity is plausible, 
as metals depleted on to dust are unlikely to have been 
returned to the gas phase within the '--^100,000-yr life- 
time of the remnant. If the metallicity is '^1/3 solar, 
this scenario would favor z 1 kpc for the arc, as higher 
remnants would be too faint and too small. 

The radial width of the arc is underpredicted by our 
SNR model. The Vela SNR, a known remnant in the 
Galactic disk, appears to be at a similar stage in its evolu- 
tion as our arc SNR model. Like the arc, its X-ray-bright 
rim is much wider than that predicted by our simulations. 
We suggest that this discrepancy could be due to the 1-D 
nature of the simulations, as well as assumptions of uni- 
form ambient density and non-thermal pressure. Higher- 
dimensional hydrodynamical simulations may better ex- 
plain the morphology of the arc, but are beyond the scope 
of this study. 

It should be noted that we have not tested other sce- 
narios for the arc's formation, and also that equilibrium 
models provide good fits to t he Suzaku spec tra. How- 
ever, we conclude by echoing fSheltonl (|2006l ). and note 
that the extraplanar SNR scenario discussed here pro- 
vides a good explanation for an arc-shaped enhancement 
in the high-latitude soft X-ray background. If this sce- 
nario is correct, it supports the idea that extraplanar 
supernovae contribute to the heating of the hot halo gas. 



We would like to thank Steve Snowden for help- 
ful comments regarding HE AS ARC 's X-ray Background 
Tool, Liz Raley for helpful discussions on fitting the 



16 



HENLEY AND SHELTON 



SNR spectral models to the data and for supplying 
the 3-D SNR simulation mentioned in !j5221 and Kyu- 
jin Kwak for calculating X-ray spectra from said sim- 
ulation. We also thank the referee, whose comments 
have helped improve this paper. This research has 
made use of data obtained from the Suzaku satellite, 
a collaborative mission between the space agencies of 



Japan (JAXA) and the USA (NASA). The OMNI data 
were obtained fro m the G SFC/SPDF OMNIWeb in- 
terface at http://omniweb.gsfc.nasa.gov. This research 
was funded by NASA grants NNH06ZDA001B-SUZ2 and 
NNX08AZ83G, awarded through the Suzaku Guest Ob- 
server Program. 



REFERENCES 



Anders, E., &; Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101: Astronomical Data 
Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 
17 

Aschenbach, B., Egger, R., & Triimper, J. 1995, Nature, 373, 587 
Balucinska-Church, M., & McCammon, D. 1992, ApJ, 400, 699 
Breitschwerdt, D., & Schmutzler, T. 1994, Nature, 371, 774 
Burrows, D. N., & Mendenhall, J. A. 1991, Nature, 351, 629 
Chen, L.-W., Fabian, A. C, & Gendreau, K. C. 1997, MNRAS, 
285, 449 

Cordes, J. M., Weisberg, J. M., Frail, D. A., Spangler, S. A., & 

Ryan, M. 1991, Nature, 354, 121 
Cravens, T. E. 2000, ApJ, 532, L153 

Cravens, T. E., Robertson, I. P., & Snowden, S. L. 2001, JGR, 106 
(All), 24883 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 
Ferriere, K. 1998, ApJ, 497, 759 
Finkbeiner, D. P. 2003, ApJS, 146, 407 
Fujimoto, R., et al. 2007, PASJ, 59, S133 

Galeazzi, M., Gupta, A., Covey, K., & Ursino, E. 2007, ApJ, 658, 
1081 

Gaustad, J. E., McCuUough, P. R., Rosing, W., & Van Buren, D. 

2001, PASP, 113, 1326 
Henley D. B., & Shelton, R. L. 2008, ApJ, 676, 335 
Ishisaki, Y., et al. 2007, PASJ, 59, S113 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 

Bajaja, E., Morras, R., & Poppel, W. G. L. 2005, A&A, 440, 775 
Koutroumpa, D., Acero, F., Lallement, R., Ballet, J., &: 

Kharchenko, V. 2007, A&A, 475, 901 
Koutroumpa, D., Lallement, R., Kharchenko, V., Dalgarno, A., 

Pepino, R., Izmodenov, V., & Quemerais, E. 2006, A&A, 460, 

289 

Koyama, K., et al. 2007, PASJ, 59, S23 

Kuntz, K. D., & Snowden, S. L. 2000, ApJ, 543, 195 

Lei, S., Shelton, R. L., & Henley, D. B. 2009, ApJ, in press 

(arXiv:0906.1532vl) 
Mitsuda, K., et al. 2007, PASJ, 59, SI 
Norman, C. A., & Ikeuchi, S. 1989, ApJ, 345, 372 
Raley E. A., Shelton, R. L., & Plewa, T. 2007, ApJ, 661, 222 
Rasmussen, J., Sommer-Larsen, J., Pedersen, K., Toft, S., Benson, 

A., Bower, R. G., & Grove, L. F. 2009, ApJ, 697, 79 
Raymond, J. C, & Smith, B. W. 1977, ApJS, 35, 419 
Robertson, I. P., & Cravens, T. E. 2003a, JGR, 108 (AlO), 8031 



— . 2003b, GeoRL, 30(8), 1439 

Savage, B. D., & Sembach, K. R. 1996, ARA&A, 34, 279 
Savage, B. D., et al. 2003, ApJS, 146, 125 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Scab, C. G. 1987, in Interstellar Processes, ed. D. J. HoUenbach & 

H. A. Thronson, Jr. (Dordrecht: Reidel), 491 
Sembach, K. R., & Savage, B. D. 1992, ApJS, 83, 147 
Shapiro, P. R., & Field, G. B. 1976, ApJ, 205, 762 
Shelton, R. L. 1998, ApJ, 504, 785 
— . 1999, ApJ, 521, 217 
— . 2006, ApJ, 638, 206 

Shelton, R. L., Sallmen, S. M., & Jenkins, E. B. 2007, ApJ, 659, 
365 

Smith, R. K., et al. 2007, PASJ, 59, S141 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 

2001, ApJ, 556, L91 
Snowden, S. L., CoUier, M. R., & Kuntz, K. D. 2004, ApJ, 610, 

1182 

Snowden, S. L., Egger, R., Finkbeiner, D. P., Freyberg, M. J., & 

Plucinsky P. P. 1998, ApJ, 493, 715 
Snowden, S. L., et al. 1997, ApJ, 485, 125 

Snowden, S. L., Freyberg, M. J., Kuntz, K. D., & Sanders, W. T. 

2000, ApJS, 128, 171 
Snowden, S. L., et al. 1995, ApJ, 454, 643 

Snowden, S. L., Mebold, U., Hirth, W., Herbstmeier, U., & Schmitt, 

J. H. M. M. 1991, Science, 252, 1529 
Spitzer, L. 1978, Physical Process in the Interstellar Medium (New 

York: Wiley) 

Taylor, J. H., Manchester, R. N., & Lyne, A. G. 1993, ApJS, 88, 
529 

Toft, S., Rasmussen, J., Sommer-Larsen, J., & Pedersen, K. 2002, 

MNRAS, 335, 799 
Winkler, P. F., Gupta, G., & Long, K. S. 2003, ApJ, 585, 324 
Yan, M., Sadeghpour, H. R., & Dalgarno, A. 1998, ApJ, 496, 1044 
Yao, Y., Nowak, M. A., Wang, Q. D., Schulz, N. S., & Canizares, 

C. R. 2008, ApJ, 672, L21 
Yao, Y., & Wang, Q. D. 2005, ApJ, 624, 751 
— . 2006, ApJ, 641, 930 
— . 2007a, ApJ, 658, 1088 
— . 2007b, ApJ, 666, 242 

Yoshino, T., et al. 2009, PASJ, submitted (arXiv:0903.2981vl) 



